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Abstract 

An online power management strategy (PMS) for plug-in 
hybrid electric vehicles is presented in this paper. This PMS 
with the capability to reduce both the fuel consumption and 
the emission, has simple mathematics and excludes a priori 
knowledge of the driving cycle. The only required 
information is the driving duration that can be estimated by 
the driver or by the vehicle information systems, so the 
proposed method can be easily implemented. Furthermore, 
an adaptive form of this PMS is presented and its 
performance is compared with other strategies. Using the 
online adaptive PMS method, the incoming driving cycle 
condition is predicted by the vehicle past conditions. In this 
paper, the engine fuel characteristics are linearized to several 
zones. At any instant, one of these zones is selected for the 
engine operation. In each zone, an optimal cost function is 
minimized for the fuel consumption and the emission 
reduction. Moreover, different cost functions are defined 
and used on various engines. Finally, the proposed PMS is 
simulated in the ADVISOR environment and compared with 
conventional method. 
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Introduction 

Increasing fossil fuel consumption has caused 
significant problems for both governments and 
societies. These non-renewable resources are depleting 
very fast. Furthermore, the air pollution and the 
greenhouse effects of their emissions, damage the 
general health. Since a big share of the fuel 
consumption is used by the conventional vehicles, the 
hybrid electric vehicles (HEVs) are suggested as an 
effective solution to reduce the fuel consumption. In 
HEVs, an electric power source is integrated into a 
conventional internal combustion engine which is 
downsized. However, the main source of their traction 
power is still the gasoline engine. The hybridization 
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and downsizing allow for fuel economy enhancements 
(Katrasnik et al., 2007). Thanks to the presence of an 
electric power train, these vehicles have the 
regenerative braking ability that increases the vehicle 
efficiency (Yu-shan Li et al., 2009). 

The latest generation of the HEVs is the Plug-in hybrid 
electric vehicles (PHEVs) that employ the grid electric 
energy. For example, they can be plugged into a 
power grid or a residential photovoltaic system 
(Gurkaynak et al, 2009). The PHEVs consume 40% to 
80% fuel than conventional HEVs (Fajri et al., 2008). 

IEEE-USA Energy Policy Committee defines that a 
PHEV should have at least a 4 kWh battery storage 
system and ten miles only electric mode driving 
distance (Shams-Zahraei et al., 2009). Therefore, the 
PHEVs have a medium storage system that is charged 
externally. Usually, the electric energy is cheaper than 
the gasoline (four times in USA (Romm et al., 2006)). 
Therefore, it is essential to use the stored energy before 
the end of the trip that makes the PMS more 
complicated. The PMS defines the component's power 
share. Many PMSs are proposed and classified into the 
HEVs and the PHEVs. From the mathematical point of 
view, the PMSs are classified into rule based and 
optimization based methods (Bayindir et al. 2011; 
Wirasingha et al. 2011). 

There are many heuristic and fuzzy rule based PMSs 
in recent literatures (Gao et al., 2010; Hui et al., 2008; 
Banvait et al., 2009; Mapelli et al., 2009; Li et al., 2010; 
Chen Zheng et al., 2009; Xiong et al. 2009). A novel 
rule-based PMS for the PHEVs that focuses on all 
electric range and charge depletion range operations is 
presented in (Gao et al., 2010). An engine on-off rule 
based control strategy with consideration on position 
of acceleration pedal is proposed in (Hui et al., 2008). 
A heuristic solution to a parallel and series-parallel 
PHEV is proposed in (Banvait et al., 2009). In this 
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paper, the energy management optimizes engine 
operational efficiency while maintaining battery state 
of charge. A rule based fuzzy logic control strategy for 
a parallel hybrid electric city public bus is proposed in 
(Li et al., 2010). A series-parallel structure for a hybrid 
electric bus is presented in (Xiong et al. 2009). In this 
paper, a fuzzy logic PMS, which consists of two fuzzy 
modules used to determine the operation mode and 
distribute torques at parallel mode respectively, is 
studied by numerical simulation. However, the 
emission reduction is exclusive of consideration in 
(Xiong et al. 2009). The rule based PMSs are simple, 
applicable, and can be calculated easily. However, the 
optimization based PMSs are usually more efficient 
than the rule based PMSs (Bayindir et al., 2011; 
Wirasingha et al. 2011). 

Generally, the optimization based PMSs are more 
complicated and accurate than the rule based PMSs 
and for both, a priori knowledge about the driving 
cycle is essential. This essential priori knowledge can 
be delivered by the solutions provided in literatures. 
However, they can lead to some more additional 
computational burden. Therefore, these PMSs are 
unsuitable for the real applications. The optimization 
methods define and minimize a cost function 
(Wirasingha et al. 2011). Both of the fuel consumption 
and the emission production are minimized easily by 
means of the proper cost function definition. The 
optimization group includes a wide spectrum of 
different methods such as the static optimization, 
numerical optimization, equivalent consumption 
minimization strategy (ECMS) and analytical 
optimization methods (Sciarretta et al., 2007). 

One of the interesting optimization PMS is dynamic 
programming (DP) (Shen et al., 2010; Yan et al., 2010). 
DP is significant time consuming with heavy 
mathematical burden. Therefore, a two-scale DP is 
proposed to solve the mentioned problems (Gong et 
al., 2007; Gong et al. 2008; Gong et al., 2009). 
Furthermore, the driving cycle prediction can be done 
by the new geographical systems like GPS, GIS, and 
traffic information system (Zhang et al., 2010; Abdul- 
Hak et al., 2009; Gong et al., 2008; Gong et al., 2007). 
(Moura et al., 2011) uses a stochastic dynamic 
programming to optimize PHEV power management 
over a distribution of drive cycles, rather than a single 
cycle by means of Markov chains. An optimization 
PMS by using the Pontryagin's minimum principle for 
the PHEVs is proposed in (Stockar et al., 2010). The 
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equivalent fuel consumption strategies (ECMS) are 
proposed by (Mapelli et al., 2009; Tulpule et al., 2009; 
Zhang et al., 2010). An artificial neural networks 
monitoring and fault diagnosis system for electric 
vehicles is proposed by (Kalogirou et al., 2000) that can 
be improved for energy management optimization. 

This paper presents a PMS that is inspired from 
(Kessels et al., 2008; Koot et al, 2005), and ECMS 
method. The presented PMS is designed for the fuel 
consumption and the emission reduction. The 
different forms of this PMS and their effects are 
discussed in this paper. 

The proposed PMS is an applicable and accurate 
method with the simple mathematic equations. 
Therefore, it has the advantages of the optimization 
and rule based methods and can be done online. This 
strategy just needs the driving duration that can be 
easily estimated by the driver or the geographical and 
information systems. 

This paper is organized as follows: An adequate 
vehicle model is presented in Section 2. The 
background of the proposed PMS and its mathematics 
are discussed in Section 3. The zone selection 
procedure is discussed in Section 4. In Section 5, an 
adaptive form of the proposed control strategy is 
discussed, and the simulation results are presented in 
Section 6. Finally, conclusions are stated in Section 7. 

Vehicle Model 

The focus of this paper is on the parallel PHEV 
topology. The schematic of this topology is shown in 
Fig. 1. In this research, the dynamics of the systems 
with low frequency (lower than one second) are 
ignored as they have insignificant effects on the fuel 
consumption of the PHEV. The proposed PMS is 
simulated in the ADVISOR environment, in which this 
tool is introduced in (Markel et al., 2002). However, 
the basic equations of this model are discussed as 
follows. 

Internal Combustion Engine (ICE) 

The fuel consumption and emission rates of an ICE are 
functions of the engine output power (Pm) and the 
engine speed as (l)-(4): 

Juel mte =f(P m ,(o) (1) 
HC mte =hc(P m ,co) (2) 
CO mte =co(P m ,(o) 
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NOx=nox(P m ,a>) 



(4) 
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FIG. 1 THE PARALLEL PHEV TOPOLOGY 

Where, <y is the engine speed and Pm is the engine 
power which is a positive parameter. In this paper, HC, 
CO, and NOx are considered as engine emissions. 

Clutch (CL) 

The vehicle engine is engaged and disengaged from 
the drive train by a clutch. The clutch will be 
disengaged during the shifting time of the gearbox, 
when the engine is off, and when the speed of the 
clutch becomes lower than the engine idle speed. 

Motor/ Generator (MG) 

There is at least one electrical machine in a PHEV that 
can act as a motor or a generator. Thus, the electric 
machine is addressed in this paper as motor/generator 
(MG). 

Unlike ICE, the output power of the MG can be 
positive or negative. In this paper, the sign of Pe is 
considered positive when the MG is operated as a 
generator and negative when it is operated as a motor. 

Neglect of friction loss at no load, the MG can be 
modelled as follows: 



1 



>gm (5) 

nmm and r]gm are the efficiency of MG in motor and 
generator modes, respectively. 



I gm 



(6) 



(7) 



Pem is the MG power on the mechanical side and Pe is 
that on the electrical side. The defined MG power 
limits in (8) are applied to the mechanical side. 



(co) < P em (a>) < P e , 



(8) 



Drive Train (D) 

The traction power (Pt) is the required wheels power 
that can be expressed as (9) (Gao et al., 2010). 

P, = v(Mgf r + ^p a C D A f v 2 +MS^ + Mgi) (W) 
2 at 

where, M is the vehicle mass; v is the vehicle speed; g 
is the gravity acceleration, p a is the air mass density; 
Cd is the aerodynamic drag coefficient of the vehicle; 
Af is the front area of the vehicle; 5 is the rotational 
inertia factor; dv/dt is the acceleration, and i is the 
grade of the road. 

If Pi is regarded as the input power to the drive train, 
it can be calculated by (10). 

P d =a d xT d (10) 

where, aid and Td are the crankshaft speed and torque, 
respectively. 



(11) 



where fd is the drive train ratio, r is the wheel radius 
and rjd is the efficiency of the drive train. 

Power Split (PS) 

The power split is assumed to have no energy losses 
and provides the following power balance: 



m 



P d +R 



em 



(13) 



Battery (B) 



The battery is modelled by a voltage source and a 
series resistance that represents the battery's internal 
and terminal losses. Therefore, (14) shows the relation 
between terminal power (Pi>) and the net internal 
power (Ps) of the battery, rf and rj- are the efficiency of 
the battery during charge and discharge modes, 
respectively. 



P h =rmx(^-,r, P s ) 
1 



The battery storage energy (£ s ) is calculated as: 
E s (t) = E s (0)+\'p s (t)dt 

JO 



(14) 



(15) 



The SOE of the battery is defined as the percentage of 
the remained energy in the battery to the total 
theoretical energy capacity of that (Eea P ). 
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SOE(t) = --^xlOO [%] (16) 

^ cap 

If the battery open clamp voltage is considered almost 
constant, the battery state of charge (SOC) can be 
regared equal to SOE. 

Electric Load (EL) 

Electric load (EL) represents auxiliaries such as head 
lights, radio, etc. The power consumed by EL is called 
PL. For simplicity, this power is assumed to be 
constant and equal to 700 W. 

Problem Definition 

It is expected that PMS maintain the vehicle 
performance and at the same time reduce the fuel 
consumption and emission. The strategy should be 
practical and easily applicable on the vehicle. 
Furthermore, the PMS should be causal. However, 
some predictive control strategies have been 
developed with some priori knowledge, but it can lead 
to some additional computational burden. The 
presented PMS has a little dependence on a priori 
information. Essential priori information for the 
proposed PMS is just the driving durationthat can be 
estimated by the assumption on the constant 
acceleration and deceleration rate as well as the speed 
limits for each road (Gong et al., 2007; Gong et al., 
2008). The proposed PMS with simple mathematics 
can be supported by common processors. Therefore, 
the presented PMS is practical and can be applied 
easily to the vehicle. 

The proposed PMS is inspired from the concepts of the 
rule based, and the optimization based strategies. 
Consequently, it has advantages of both. 

This strategy is mainly based on ICE characteristics. In 
the vehicles, the ICE consumes fuel and produces 
pollution. Therefore, the proper ICE handling causes 
the significant emission and fuel consumption 
reduction. Generally, in many literatures, the engine 
torque-speed map is considered. However, this study 
is established base on the fuel consumption 
characteristics. 

The fuel consumption characteristics for SI41 and 
Prius ADVISOR engine models are shown in Fig. 2 
and Fig. 3, in which the curves are related to different 
engine speeds. Each curve is linearized and discretized 
into several zones that limit the engine operation. In 
this paper, each curve is divided into 11 zones. The 
SI41 has a very rugged fuel characteristic, whereas the 



fuel characteristic of Prius is monotonous. In other 
words, the slope of the Prius curves decreases with the 
increase of power. The Prius engine special structure 
results from this monotonous characteristic. The Prius 
engine torque-speed map is shown in Fig. 4in which 
that at a constant speed, the engine efficiency increases 
due to the rise in torque causes. However, the most 
efficient area of common engines is usually in the 
areas lower than the maximum torque. For instance, 
the torque-speed map of the SI41 engine which is a 
common engine is shown in Fig. 5. 
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FIG. 2 FUEL MAP CHARACTERISTIC FOR THE SI41 ENGINE 
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FIG. 3 FUEL MAP CHARACTERISTIC FOR THE PRIUS ENGINE 
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FIG. 5 SI41 FUEL CONSUMPTION MAP IN GRAM PER KWH 

A power management strategy for HEVs is presented 
in (Kessels et al., 2008). The idea of (Kessels et al., 2008) 
is that the engine should work in the area with high 
(or low dependence on the engine map) slope when 
the SOC of the battery is low and work in the area 
with low (or high) slope when the SOC is high. 
Therefore, PMS can be compatible for engines like the 
Prius, but it is unsuitable for engines with a rugged 
fuel map like SI41. 

Fig. 2 and Fig. 3 show that the fuel consumption 
decreases if the engine speed descends with constant 
power. Therefore, in series PHEVs or in the topologies 
that the speed of the engine is controllable; the engine 
operation at low speeds will be suitable. 

The slope of each zone is an important parameter of 
the proposed PMS. Typically, the slope of each zone 
(A//) is constant, where i indicates the zone number. A/; 
represents the extra fuel mass flow needed to produce 
a small amount of mechanical power by the engine. 



(17) 

The slopes of different zones and curves are calculated 
and gathered in a lookup table. Furthermore, the 
similar lookup tables are prepared for emission 
characteristics. The characteristics of these emissions 
for SI41engine are shown in the Figs. 6-8. 

The variations of HC and CO against power changes 
are similar to the fuel usage variations. Consequently, 
the attempt on fuel reduction causes a reduction of CO 
and HC emissions. However, the difference between 
the NOx and the fuel characteristics is significant. In 
fact, there are some zones with high NOx emission in 



the optimal fuel consumption area. This difference can 
have important effects on the PMS for NOx reduction. 
Hence, for reduction of the NOx emissions, these zones 
can be avoided. Certainly, the fuel consumption rate is 
increased, but it is not significant. In fact, the engine 
operating point is still in the optimal fuel usage area , 
inside of which there are just some holes. Thus, a 
suitable balance between fuel usage and NOx emission 
is achieved. 

Different engines have various fuel and emission 
characteristics. Consequently, this strategy 
performance and its setting will be different from each 
other. There is more discussion about the effect of 
engine characteristics on the PMS results in Section 6. 

The emission curves should be discretized and 
linearized same as fuel consumption curves. The slope 
of each zone of the HC, CO and NOx characteristic 
curves are given in (18)-(20), respectively. 
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FIG. 6 HC EMISSION MAP CHARACTERISTIC FOR SI41 ENGINE 




FIG. 7 CO EMISSION MAP CHARACTERISTIC FOR SI41 ENGINE 
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FIG. 8 NOx EMISSION MAP CHARACTERISTIC FOR SI41 
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dhc(P m ,co) 

dco(P m ,co) 
dP 

m 

dnox(P m ,co) 
dP 



(18) 



(19) 



Thus, (l)-(4) can be linearized as follows: 
f.iP^^Kf^ + A^co)^ 
he, (P m , co) * hc 0i (CO) + A hci (P m , CD).P m 
co, (P m , co) * co 0i (co) + A cm (P m , co).P m 
nox t (P m , co) « nox 0i (co) + k noxi (P m , co).P m 
i = \,2,...,N f 



(20) 

(21) 
(22) 
(23) 

(24) 

where fa, hcoi, coo, , and noxoi indicate the initial fuel 
usage and emission rates of the HC, CO, and NOx for 
the ;-th zone, respectively, and Nf is the number of the 
zones. 

By the selection of the proper zone, the approximate 
fuel consumption and emission can be calculated. In 
order to find this zone, it is required to find the 
optimal value of A/ that is explained in more details in 
Section 4. Then by solving an instantaneous optimal 
cost function, an appropriate power share value for 
the ICE and MG is calculated. This instantaneous 
optimal cost function is derived from the cost function 
of the ECMS method, in which an instantaneous cost 
function is defined for equivalent fuel minimization 



(Mapelli et al., 2009). (25) explains this cost function 
where J(t) is the equivalent fuel usage, f(Pm) is the real 
fuel usage by ICE, Hi is the low heating value of the 
fuel, and s is the equivalent factor that depends on the 
driving cycle and the efficiency of components of the 
power train (Zhang et al., 2010). 



J(t) = f(PJ- 



H, 



(25) 

By linearizing the fuel map, the ratio of the equivalent 
factor and Hi, is considered equal to the slope of the 
chosen operation zone (A/). For each zone (i), there is a 
certain slope called A/;. Moreover, in order to acquire 
only one variable, Pm can be replaced with (13). 
Therefore, (25) is rewritten as: 

J(t) = f(P em \P li ,P l ,co)-A fi P em (26) 

The cost function can be formed for the emissions. To 
achieve a thorough cost function, all of these equations 
are added together by the weight coefficients. (27) is 
the comprehensive cost function where Wf, Whc, Ww, and 
Wnax are the weight coefficients that represent the 
importance of each part of the cost function. 

J(t) = w f (f(P em \P d ,P„ ( o)-A fi PJ + 

^(hc(P em \P d ,P„co)-\ a PJ + 
wJco(P m \P d ,P l ,CD)-2, coi PJ+ (27) 
KoMox(PJP d ,P n co)-\ oxl PJ 

The proper zone selection is based on the fuel curves. 
Therefore, the engine operation in the optimal fuel 
consumption area is guaranteed based on the 
proposed PMS. In other words, if the cost function just 
covers the emissions, the fuel consumption would not 
increase significantly. Consequently, in this paper, the 
values of the weight coefficients for emission 
reduction are assumed to be wj=0, u>h<=\, wm=\ and 

Wnox=l. 

In this study, two control variables (Pem and S) are 
used to manipulate the ICE and MG. Pem is the power 
of MG in the mechanical side, and S is the binary 
variable that defines on or off mode of the ICE. If S 
equals 1, and the engine speed is more than its idle 
speed, then the ICE will be on. Otherwise, it will be off. 

S will be zero in the following situations: 

• If SOC of the battery is more than SOCl and at 
the same time, the MG and batteries can 
provide enough power for traction. In other 
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words, if (28) and (29) are satisfied. 

-P, > P 

1 a — 1 " 



em_rmn 



(28) 



P L - P b van 



7 mm (29) 
• In the regenerative braking mode, if the ICE is 
off, the kinetic energy of the wheels will be 
captured by MG. It is assumed that just 60% of 
this energy can be regenerated (Kalogirou et al., 
2008). 

The results should satisfy the following limits. 
P (t)< P (t)< P ( 

1 em_mm V 1 / — 1 em\ l > — 1 em_max 1 
P s ^(t)<P s (f)<P s 



:(t) Vfe[0,fj 

s ,nax(0 Vfe[0,f e ] 



Zone n 



(30) 
(31) 
(32) 
(33) 



Tnn(0^ Zone{t) < Zone^ (t) V? £ [0, t e ] 
SOC(t e )*SOC L 

SOC(t) > SOC L Vf e (0, t e ) , 34 s 

However, the ICE can work in higher zones compared 
to the optimal area, if the SOC is lower than SOCl or 




when the driving cycle is missing, these constraints 
can be defined by design engineers. 

Online Zone Selection 

The zone selection is an important section of the 
proposed strategy. The two variables (dSOC and A*) 
are used to manipulate the zones. dSOC indicates the 
SOC situation such as (35) where, SOCref is reference 
value for the battery SOC. The zone shifting direction 
is determined by dSOC and time of the zone shifting is 
dependent on the value of A*. 



SOC 



ref 



(35) 



In each loop of the program, the zone number is 
allowed to change just one unit. In other words, the 
engine operating point can jump from present zone 
just to the neighbor zones. For positive dSOC, it will 
jump to the zone with lower power. Otherwise, the 
zone shifting direction will be towards the neighbor 
zone that has higher power points. 
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FIG. 9 FEEDBACK DIAGRAM FOR ESTIMATING A* 



The A* is the online calculated parameter to select the 
appropriate zone for the engine operating point. The 
zone shiftingoccurs when the slope of the neighbor 
zone is closer to the A* than the present zone slope. 
The A* is calculated by a system shown in Fig. 9. There 
is a limitation for A* which depends on the 
characteristics of the engine map. The A* is a 
comparative parameter between the slope of the 
neighbor zones. Therefore, the lower and upper 
constraints of the A* are the minimum and maximum 
members of the fuel slope lookup table, respectively. 

The input signal of the proportional integrator (PI) 
controller is multiplied by -1, if the slope of the 
neighbor zones has reduced with the increasing power. 
Otherwise, the multiplier is +1. So, the A* is equal to 



A> — An + 



Kp.eit) + Kj f e(v)dv 
Jo 



(36) 

where Ao is an initial guess, e(t) is the input of the PI 



controller. Selection of the parameters Kp and Ki for a 
small closed-loop bandwidth allows tracking the 
SOCref by means of the actual SOC (Kalogirou et al., 
2008). 

The proper initial value for A* (Ao) is the slope of the 
most optimum zone. With this selection, the engine 
operates near the efficient zone at the start of the trip. 
Appropriately, tuning of the parameters Kp and Ki is 
described and analyzed in (Kessels et al., 2007). The 
considered values for Ao, Kp and Ki are 25 x 10-7, 67><10-8, 
and 33 x 10-5, respectively. 

The SOCnf has an essential role in choosing the proper 
zone. The DP SOC trajectory for the predictive driving 
cycle is assumed to be an appropriate trajectory for 
SOC, therefore, in the literatures the DP is used as a 
comparative method (Yang et al., 2010; Gong et al., 
2007; Kessels et al., 2008). However, this method is 
time consuming and has heavy numerical calculations. 
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Therefore, an easy and practical SOCref is suggested. 
The SOCnf is considered as a straight line from SOCh at 
the beginning of the trip to SOCl at the end. However, 
if the battery SOC reaches SOCl before the end of the 
trip, then the SOCref will be switched to the horizontal 
line of the SOCl value. These two situations are shown 
in Fig. 10 and Fig. 11. 

SOCref is a reference and guide to manipulate the zones 
and to calculate A*. Hence, in some low speed urban 
areas, the engine is usually off. Therefore, the MG 
provides the traction power for most of the time. 
Consequently, SOC will be lower than SOCref in the CD 
(charge depletion) mode such as shown in Fig. 10. 



0.x 



0.7 



0.6 



U 0.5 

w 

U0.4 

91 

0.3h 



0.2- 



0.1 





— SOC 
— SOCref 






\ 

- 








s 

w, 




i i i i i i 



6000 8C 
Time (sec) 



12000 14000 



FIG. 10 SOC AND SOCref FOR 10 CYCLES OF UDDS (Alipour et al., 

2012) 
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FIG. 11 SOC AND SOCref FOR 10 CYCLES OF INDIA_ HWY 
(Alipour et al., 2012) 

Fig. 11 shows the results of 10 repetitions of the 
India_HWY driving cycle, in which, speed is high; 



therefore, usually the engine is on and provides most 
of the traction power. Hence, the battery is depleted 
slowly, but it should be depleted to SOCl before the 
end of the trip. The SOCref is used for this goal and the 
length of the trip defines the slope of the SOCref. A 
correct trip length prediction has an important role in 
the battery discharging slope. 

If the driver changes his destination in the middle of 
the trip, a new SOCref is established which is a straight 
line from the present SOC to SOCl at the end of the 
trip. If the driver does not specify the destination, the 
SOCnf can be a horizontal line with SOCl value. 
Therefore, the engine will be off in charge depletion 
mode. However, this condition is not optimal. 

Online Adaptive Strategy 

The proposed PMS fails to consider the driving cycle 
condition. However, the vehicle efficiency may be 
increased duo to the consideration of the driving 
pattern. In this paper, the driving cycles are classified 
into two groups as urban and highway. Generally, in 
the highway, the power demand and the vehicle speed 
are high and the engine can easily work in the efficient 
area. However, in the urban areas the vehicle speed is 
low, and it stops working frequently. Furthermore, it 
is possible that a driver drives his vehicle in different 
driving patterns in one day. Therefore, it is logical that 
the battery can be charged in the highway and while 
discharged in the urban. 

This idea is added to the proposed PMS by shifting the 
operating zone to the neighbor zone with higher 
power and keeping the engine on (if SOC is lower than 
SOCref) in the highway driving cycles. In the urban 
cycles, the engine should be operated in the low power 
zones. 

The driver can manually define the type of the driving 
conditions (highway or urban). However, each cycle 
can be divided into several small highways and urban 
cycles. It is assumed that the road condition ahead will 
be similar to the recent past condition (say 40 seconds 
for these durations). So, by checking the driving cycle 
of the past moments, the road condition is defined in 
the future. If the vehicle speed at all moments of a past 
period or the average speed of the vehicle is greater 
than a certain speed value, then the next period should 
be considered as a highway, otherwise it is considered 
as urban. This specific speed is acquired by trial and 
error. In this paper, the used specific speed is 32.4 km/h. 
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In this strategy, dsoc is the adapted factor. Therefore, 
the driver cycle mode effects the zone selection. 
Moreover, the battery should be charged in the 
highway, and be used in the urban. This online 
adaptive strategy can improve the fuel saving mainly 
in the combined driving cycles, but it may increase the 
emissions. In the highway driving cycles, the engine 

Time 



operating points are moved rapidly towards the high 
power areas. Generally in these areas, the emissions 
are not in the optimum points. However, this subject is 
dependent on the engine characteristics, but it is the 
case in most of them. 

The implementation schematic for the proposed PMS 
is shown in Fig. 12. 
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FIG. 12 THE IMPLEMENTATION SCHEMATIC FOR THE PROPOSED PMS 

TABLE 1 THE SIMULATED VEHICLE PARAMETERS 



The proposed PMS is simulated in the ADVISOR 
environment for several driving cycles (Alipour et al., 
2010). The results are compared with that of a heuristic 
rule based PMS proposed in (Banvait et al., 2009). The 
assumed vehicle parameters are presented in Table 1. 
The components of this PHEV are given in Table 2. 
Furthermore, another simulation is performed for this 
vehicle with the SI63 engine. The vehicle weight with 
this new engine is 1470 kgs. 

The simulation is performed for five combined driving 
cycles and ten repeats of seven normal driving cycles. 
A list of these cycles is given in Table 3. The five 
combined driving cycles are formed by different 
highway and urban driving cycles. The average of the 
fuel consumption and emissions for these cycles is 
considered as a comparison criterion. However, a 
discussion is performed on the simulation results for 
ten repeats of the NEDC driving cycle. 



Quantity 


Notation 


Value 


unit 


Vehicle mass 


M 


1400 


K g 


Front area 


Ad 


2.0 


m 2 


Air drag 


a 


0.3 


- 


Rolling 
resistance 


Cr 


0.015 




Air density P 1.2 Kg/m 3 


Gravity 


G 


9.8 


mis 2 


Wheel radius 




0.3 




Upper SOC level 


SOCh 


0.8 




Lower SOC level 


SOCl 






Gear ratios 


8' 


13.45_7.57_5.01_3.77_3.01 




TABLE 2 THE SIMULATED VEHICLE COMPONENTS 



Engine 



Battery 



SI41, with 41 kW maximum power and 0.34 

peak efficiency, 1 .0L, 3cylinders. 
AC59, with 56 kW maximum power and 0.91 
peak efficiency 
NIMH28_OVANIC, with (C/3) 28 Ah 
nominal capacity and 335 V nominal voltage, 

6V cell nominal voltage, 3.6kg weight, 175 
Wh nominal energy (C/3), 1.6kW peak power 
(10s pulse) 
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TABLE 3 LIST OF THE DRIVING CYCLES 



Table 4. 



Seven normal driving 
cycles 


Five combined driving cycles 


Cyc_EUDC 
Cyc.FTP 

Cyc_HWFET 
Cyc_India_urban_sample 
Cyc_India_hwy_sample 

CycJMEDC 

CycJJDDS 


(India)_ 
hwy_urban_hwy_urban 
(India)_urban_hwy_hyw_urban 

(India)_urban_hwy_urban 
(India)_hwy_hwy_urban_urban 
UDDS_HWFET_UDDS (2cycle) 



The simulations are performed for three different 
strategieswhich are the proposed PMS with and 
without adaptation, and the rule based strategy from 
(Banvait et al., 2009). The expected goals of these 
strategies are the fuel, fuel and emission, emission, 
and NOx reduction. Furthermore, the effects of the 
zone correction for the NOx reduction that has been 
introduced in Section 5 are studied. The simulation 
results for seven normal driving cycles are calculated. 
Then the average value for these cycles is computed 
for each condition. The average values are shown in 

TABLE 4 THE AVERAGE VALUES OF THE SIMULATION RESULTS FOR THE SI41 



The results for the non-adaptive PMS with the fuel 
optimization cost function are considered as base 
values for comparison. The fuel consumption is 
increased by changing the cost function from the fuel 
reduction to the emission reduction. 

The SI41 NOx emission characteristic is very different 
from the fuel usage. Consequently, the most increase 
in the fuel consumption is in the NOx reduction case. 
Particularly, the fuel consumption is increased when 
the engine operation area is corrected by taking into 
account the NOx map. However, the NOx emission 
will be minimized in this case. There is 2.8% less NOx 
emission for the non-adaptive PMS with the zone 
correction in comparison to the PMS without the zone 
correction. At the same time, it has 0.6% more fuel 
consumption. There are similar results for the adaptive 
PMS. Consequently, the zone correction method is 
more beneficial. Specially, with considering that the 
NOx reduction is the goal of the optimization. 

ENGINE (AVERAGE PERFORMANCES OF THE SEVEN NORMAL CYCLES) 



Power Management Strategy 


Fuel Usage 


HC 

Emission 


CO 

Emission 


NOx 

Emission 


Total 
Emission 


Liter/lOOKm 


% 


g/s 


% 


g/s 


% 


g/s 


% 


g/s 


% 


Proposed PMS without adaption-with the cost 
function for the fuel reduction 


3.57 


100 


0.100 


100 


0.275 


100 


0.141 


100 


0.516 


100 


Proposed PMS without adaption-with the cost 
function for the emission reduction 


3.64 


101.9 


0.098 


98 


0.271 


98.5 


0.126 


89.3 


0.495 


95.9 


Proposed PMS without adaption-with the cost 
function for fuel and emission reduction 


3.64 


101.9 


0.098 


98 


0.271 


98.5 


0.126 


89.3 


0.495 


95.9 


Proposed PMS without adaption-with the cost 
function for NOx reduction with zone correction 


3.67 


102.8 


0.098 


98 


0.294 


106.9 


0.121 


85.8 


0.513 


99.4 


Proposed PMS without adaption-with the cost 
function for NOx reduction without zone correction 


3.65 


102.2 


0.098 


98 


0.278 


101.1 


0.125 


88.6 


0.501 


97.1 


Proposed PMS with adaption-with the cost function 
for the fuel reduction 


3.55 


99.4 


0.113 


113 


0.333 


121 


0.152 


107.8 


0.598 


115.9 


Proposed PMS with adaption-with the cost function 
for the emission reduction 


3.57 


100 


0.111 


111 


0.305 


110.9 


0.140 


99.3 


0.556 


107.7 


Proposed PMS with adaption-with the cost function 
for fuel and emission reduction 


3.57 


100 


0.111 


111 


0.305 


110.9 


0.140 


99.3 


0.556 


107.7 


Proposed PMS with adaption-with the cost function 
for NOx reduction with zone correction 


3.64 


101.9 


0.109 


109 


0.336 


122.2 


0.127 


90.1 


0.572 


110.8 


Proposed PMS with adaption-with the cost function 
for NOx reduction without zone correction 


3.61 


101.1 


0.111 


111 


0.337 


122.5 


0.141 


100 


0.589 


114.1 


rule based PMS from (Banvait et al., 2009) 


3.70 


103.6 


0.096 


96 


0.649 


236 


0.142 


100.7 


0.887 


171.9 
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TABLE 5 THE AVERAGE VALUES OF THE SIMULATION RESULTS FOR THE SI41 ENGINE (AVERAGE PERFORMANCES OF THE FIVE COMBINED CYCLES) 



Power Management Strategy 


Fuel Usage 


HC 

Emission 


CO 

Emission 


NOx 

Emission 


Total 
Emission 


T iter/1 f)f)Km 


% 


g/s 


% 


g/s 


% 


g/s 


% 


g/s 


% 


Proposed PMS without adaption-with the cost 
function for the fuel reduction 


2.16 


100 


0.109 


100 


0.344 


100 


0.124 


100 


0.577 


100 


Proposed PMS without adaption-with the cost 
function for the emission reduction 


Z.16 


1 nn 
1UU 


U.lUo 


96.3 


U.3U3 


QQ 1 
OO.i 


U.lUo 


O/.l 


(1 CI/ 

U.516 


QQ A 


Proposed PMS without adaption-with the cost 
function for NO* reduction with zone correction 


2.16 


100 


0.105 


96.3 


0.339 


98.5 


0.103 


83.1 


0.547 


94.8 


Proposed PMS with adaption-with the cost 
function for the fuel reduction 


2.08 


96.3 


0.144 


132.1 


0.469 


136.3 


0.170 


137.1 


0.783 


135.7 


Proposed PMS with adaption-with the cost 
function for the emission reduction 


2.12 


98.1 


0.141 


129.3 


0.444 


129.1 


0.162 


130.6 


0.747 


129.5 


Proposed PMS with adaption-with the cost 
function for NO* reduction with zone correction 


2.14 


99.1 


0.138 


126.6 


0.448 


130.2 


0.139 


112.1 


0.725 


125.6 


rule based PMS from (Banvait et al, 2009) 


2.18 


100.9 


0.106 


97.2 


0.494 


143.6 


0.135 


108.9 


0.7368 


127.4 



TABLE 6 THE AVERAGE VALUES OF THE SIMULATION RESULTS FOR THE SI63 ENGINE (AVERAGE PERFORMANCES OF THE FIVE COMBINED CYCLES) 



Power Management Strategy 


Fuel Usage 


HC 

Emission 


CO 

Emission 


NOx 

Emission 


Total Emission 


Liter/lOOKm 


% 


g/s 


% 


g/s 


% 


g/s 


% 


g/s 


% 


Proposed PMS without adaption-with the 
cost function for the fuel reduction 


2.36 


100 


0.083 


100 


0.374 


100 


0.077 


100 


0.534 


100 


Proposed PMS without adaption-with the 
cost function for the emission reduction 


2.34 


99.2 


0.083 


100 


0.372 


99.5 


0.076 


98.7 


0.531 


99.4 


Proposed PMS with adaption-with the cost 
function for the fuel reduction 


2.28 


96.6 


0.102 


122.9 


0.495 


132.4 


0.103 


133.8 


0.7 


131.1 


Proposed PMS with adaption-with the cost 
function for the emission reduction 


2.28 


96.6 


0.103 


124.1 


0.494 


132.1 


0.103 


133.8 


0.7 


131.1 



Table 4 shows that the adaptive strategy causes this strategy would have more advantage for the 
insignificant improvement in the fuel saving for the combined driving cycles. The adaptive strategy has 
normal driving cycles. However, it is expected that considerable increase in the total vehicle emission. The 
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engine operating points for the adaptive PMS move 
very fast towards higher power zones in the highway 
periods, but these areas are improper in the emissions 
point of view. 

The strategy to reduce the emission has the same 
results as the strategy with both emission and fuel 
optimization method. However, it has a simple cost 
function. Therefore, the emission cost function is a 
proper function to optimize the fuel and emission. 

In comparison with the rule based strategy, the 
presented PMS has more advantages in both the 
emission reduction, and the fuel saving. The rule 
based PMS has considerable CO emission. In this 
strategy, the engine is operated in the high torque 
areas as Fig. 13. The CO emission rate is high in this 
district for the SI41 engine. 




TOO 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 



Speed (ipm) 

FIG. 13 ENGINE OPERATING POINTS FOR THE RULE BASED 
PMS IN 10 NEDC (G/KWH) 

The maps for the fuel usage and NOx reduction PMSs 
are shown in Fig. 14 and Fig. 15, respectively. In Fig. 
15, the points inside the dashed ellipses have caused 
the high fuel consumption for the NOx reduction 
mode. In fact, the defined cost function for NOx 
reduction strategy moves the operating points towards 
the low torque areas, in which the NOx emission rate is 
low. The NOx maps for these two situations are shown 
in Fig. 16 and Fig. 17. The engine operating points in 
Fig. 17 are closer to the low NOx emission districts 
than that in Fig. 16. 

The mentioned methods are simulated for the 
combined driving cycles, for which the simulation 
results are calculated. and then as well as the average 
value for each condition. The average values are 
shown in Table 5. The fuel consumption is reduced for 
the adaptive PMS. However, the other results are 
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similar to that of the normal driving cycles. 




Speed (rpm) 



FIG. 14 THE SI41 FUEL MAP FOR THE FUEL REDUCTON COST 
FUNCTION IN NEDC CYCLE (G/KWH) 
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FIG. 15 THE SI41 FUEL MAP FOR THE NOx REDUCTION COST 
FUNCTION FN NEDC CYCLE (G/KWH) 
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500 1000 151 2000 251 3000 3500 4000 4500 501 5500 WOO 
Speed (rpm) 

FIG. 17 THE SI41 NOx MAP FOR THE NOx REDUCTION COST 
FUNCTION IN NEDC CYCLE (GIKWH) 
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FIG. 18 NOx EMISSION MAP CHARACTERISTIC FOR THE SI63 
ENGINE 

Another simulation is performed for the SI63 engine. 
The features of this engine have some differences with 
the SI41 engine. For instance, the SI63 emission and 
fuel characteristics are almost similar together. 
Therefore, the zone correction for the NOx emission 
cannot be performed for this engine. The SI63 NOx 
characteristic is shown in Fig. 18. The simulation 
results for five combined driving cycles are calculated. 
Then the average value for these cycles is calculated 
for each condition, which is shown in Table 6. There is 
a considerable increase in the fuel saving and 
emissions for the adaptive strategy, which is in 
accordance with the results of the SI41 engine. 
However, the change in the cost function has a little 
change in the fuel and emission results of the SI63 
engine whose characteristics are very important in the 
presented PMS method. 



Conclusions 

This paper has presented a powerful online power 
management strategy for the parallel PHEVs. This 
PMS excludes from the need of complex mathematics 
as well as a priori driving cycle information. Due to 
these features, the proposed PMS is an applicable 
strategyand can reduce the emissions and the fuel 
consumption as well. 

Furthermore, an online adaptive version of this PMS 
has been presented. In particular, the adaptive PMS 
has reduced the fuel consumption in the combined 
driving cycles. However, the emission has been 
increased in this case. Therefore, the use of the 
adaptive PMS is dependent on the environmental laws. 

The engine characteristics have an important role in 
the performance of this PMS. In fact, this PMS has 
been established based on the engine characteristics. 
Therefore, the emission and the fuel consumption 
reduction have varying intensity for different engines. 
However, these effects are almost similar. 
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